Quantum Monte Carlo calculations of magnetic moments and Ml 
transitions in A < 7 nuclei including meson-exchange currents 

L. E. Marcucci^'^jll Muslema Pervin^jll Steven C. Pieper^jil 
R. Schiavilla^'^i and R. B. Wiringafl 

^Department of Physics "Enrico Fermi", University of Pisa, 1-56127 Pisa, Italy 
^INFN, Seztone di Pisa, 1-56127 Pisa, Italy 
^Physics Division, Argonne National Laboratory, Argonne, Illinois 60439 
^Theory Center, Jefferson Laboratory, Newport News, Virginia 23606 
^Department of Physics, Old Dominion University, Norfolk, Virginia 23529 

(Dated: October 2, 2008) 

Abstract 

Green's function Monte Carlo calculations of magnetic moments and Ml transitions including 
two-body meson-exchange current (MEC) contributions are reported for A <7 nuclei. The realistic 
Argonne fig two-nucleon and Illinois-2 three-nucleon potentials are used to generate the nuclear 
wave functions. The two-body meson-exchange operators are constructed to satisfy the continuity 
equation with the Argonne vis potential. The MEC contributions increase the ^4=3, 7 isovector 
magnetic moments by 16% and the A=6,7 Ml transition rates by 17-34%, bringing them into very 
good agreement with the experimental data. 

PACS numbers: 21.10.-k, 23.20.-g, 23.40.-s 



Electronic address: marcucci@df.unipi.it 



^Electronic address: |muslema@phy.anl.gov| 
■l-Electronic address: |spieper @anl.gov| 
^Electronic address: schiavil@jlab.org| 
^Electronic address: 



wiringa@anl.gov 



1 



I. INTRODUCTION 



In a recent paper we reported quantum Monte Carlo (QMC) calculations of electroweak 
transitions in A = 6, 7 nuclei. The QMC method is a two-step process, with an initial 
variational Monte Carlo (VMC) calculation to find a good trial function, followed by a 
Green's function Monte Carlo (GFMC) calculation to refine the solution. When used with 
the Argonne fig two-nucleon j2] {NN) and Illinois-2 three-nucleon (3A^) potentials, the 
final GFMC results reproduce the ground- and excited-state energies for A < 10 nuclei 
0, S, B 9 very well. 

In Ref. [H we studied magnetic dipole (Ml) and electric quadrupole {E2) transitions and 
nuclear beta-decay (Fermi and Gamow- Teller) rates. These were the first off-diagonal matrix 
element calculations using the nuclear GFMC method. However, only one-body transition 
operators were used to calculate the matrix elements. We noted that two-body meson- 
exchange- current (MEG) operators are known to increase isovector magnetic moments by 
15-20% for A=3 nuclei |8|, while a previous VMC calculation for the width of the first Ml 
transition in ^Li was also increased by 20% [91]. In this paper we use GFMC wave functions 
to investigate MEG contributions to magnetic moments for the ground states of A=2-7 
nuclei as well as a number of Ml transitions in ^4=6, 7 nuclei. We find significant isovector 
contributions from the MEG operators and overall very good agreement with experiment. 

A brief review of the QMC calculational method is given in Sec. [Tll The electromagnetic 
current operator is discussed in detail in Sec. Illli Results and conclusions are given in 
Sees. IV] and El 



II. QUANTUM MONTE CARLO METHOD FOR TRANSITIONS 

We evaluate the diagonal magnetic moment matrix element {'^{J'^;T)\0\'${J'"]T)) and 
the off-diagonal transition matrix element {"^ f{J'" ]T')\0\'^i{J'"]T)) , where O is the full 
electromagnetic operator. The nuclear wave function with a specific spin-parity J'^ and 
isospin T is denoted as \E'(J''; T) and is a solution of the many-body Schrodinger equation 

Hm{r]T) = Em{r]T) . (i) 

The Hamiltonian used here has the form 



H = Y.K^ + Y.^^,+ E (2) 

i i<j i<j<k 

where Kiis the nonrelativistic kinetic energy and Vij and Vijk are respectively the Argonne fig 
(AV18) 0] and Illinois-2 (IL2) Q potentials. The VMC trial function ^t(^^;T) for a given 
nucleus is constructed from products of two- and three-body correlation operators acting on 
an antisymmetric single-particle state of the appropriate quantum numbers. The correlation 
operators are designed to refiect the infiuence of the interactions at short distances, while 
appropriate boundary conditions are imposed at long range 0, [il|. The ^t(J^;T) has 



embedded variational parameters that are adjusted to minimize the expectation value 
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which is evaluated by Metropohs Monte Carlo integration 12(]. Here Eq is the exact lowest 



eigenvalue of H for the specified quantum numbers. A good variational trial function has 
the form 



1+ E 

i<j<k 



fjTNI 



(4) 



The Jastrow wave function, is fully antisymmetric and has the {J'";T) quantum numbers 
of the state of interest, while Uij and Ujjk^ are the two- and three-body correlation operators. 
More details may be found in Ref. |lt]. The error in the variational energy Ey is of order 



1^1 



P, where \l/o is the exact lowest-energy eigenstate of H for a given set 
of quantum numbers. Other expectation values calculated with have errors of order 
|\E'o - ^t|/|^t|- 



The GFMC method [13|, [ij] reduces the VMC errors by using the relation 

^0 = lim exp[-(i7 - £'o)r]^T ; 



(5) 



that is the operator exp[— (if — -E'o)''"] projects \l/o out of "^t- If the maximum r actually 
used is large enough, the eigenvalue Eq is calculated exactly while other expectation values 



are generally calculated neglecting terms of order |\l/o ~ ^Tp/|^rP and higher (ill ]. 

In the following we present a brief overview of the nuclear GFMC method; much more 



detail may be found in Refs. [Ill 
propagated wave function \l/(r) 



15j. We start with the of Eq. and define the 



-{H~Eo)Ar 



(6) 



where we have introduced a small time step, r = riAr; obviously \&(r = 0) = and 
\l/(r — >■ oo) = \&o- Quantities of interest are evaluated in terms of a "mixed" expectation 
value between "ifx and \l/(r): 



(0(r)) 



M 



(^(r)|0|^r) 



(7) 



(^(r)I^T) 

The desired expectation values would, of course, have ^^(t) on both sides; by writing \E'(t) = 
(5\l'(r) and neglecting terms of order [(5\l/(r)]^, we obtain the approximate expression 

{^ir)\0\^{r)) 



(0(r)) 



(0(r))M + [(0(r))M - {0)v] 



(8) 



(^(r)l^(r)) 

where (0)v is the variational expectation value. 

For off-diagonal matrix elements relevant to this work the generalized mixed estimate is 
given by the expression 

($^(r)|0|^^(r)) 



{0{r))M. + {0{r))M, - {0)v 



where 



(0(r))M, 



{¥t\0\^'{t)) 



(9) 

(10) 
(11) 



and {0{T))Mf is defined similarly. For more details see Eqs. (19-24) and the accompanying 
discussions in Ref. Ili. 
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III. THE ELECTROMAGNETIC CURRENT OPERATOR 



The model used for the nuclear electromagnetic current operator j{q) is based on the 
study of Ref. [l6|. It represents j{q) as a sum of one-, two- and three-body terms that 
operate on the nucleon degrees of freedom, 



i i<j i<j<k 



(12) 



q being the three- momentum transfer. The one-body operator ji{q) is derived from the 
nonrelativistic reduction of the covariant single-nucleon current, by expanding in inverse 
powers of the nucleon mass m. In the notation of Ref. [8|, it is written as 



j.(9) = |^te,e--} + ^/i.^.xge 



iqri 



where {■ ■ 



■} denotes the anticommutator, the quantities and /ij are defined as 

1 r 



(13) 

(14) 
(15) 



and p, (T and r are the nucleon's momentum, Pauli spin and isospin operators, respectively. 
Finally, G^{qfJ {G^{qf^)) and G^{qf^) {G^{qf^)) are the isoscalar and isovector combinations 
of the nucleon electric (magnetic) Sachs form factors, respectively, evaluated at the four- 



momentum transfer qi = q^ — cu^ with 



q^ + M| — M(, where Mj and Mf are initial 



and final nuclear masses (only elastic scattering or inelastic scattering to discrete final states 
are considered in the present work). 

The current operator satisfies the current conservation relation (CCR) 



q-j{q) = [H,p{q)] 



(16) 



Here H is the nuclear Hamiltonian consisting of two- and three-nucleon interactions, the 
AV18 [2] and IL2 [3| potentials, respectively, and p{q) is the charge operator which, to 
lowest order in 1/m, is written as 



with 



pAq) 



e,;e 



(17) 
(18) 



To this order, the CCR separates into 



q-3i[q) 
Q-jij{q) 



—,piiq) 
Vij,ptiq) - 



Piiq)] 



(19) 
(20) 



and similarly for the three-body current jijk{q)- The one-body current of Eq. f|T3l) is easily 
seen to satisfy Eq. (fT9|l . 
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A. Two-nucleon current 



The two-body current operator jij{q) is separated into two parts, labeled model- 
independent (MI) and model-dependent (MD), following the scheme of Ref. [13] • The MI 
two-body currents have longitudinal components that satisfy the CCR of Eq. fl2U]) with the 
AW potential Vij, i.e. the AV18 |2|. The potential can be written as 

% = + < ' = 4 + ' (21) 

where ffj" and v^^^ are the isospin-symmetry conserving (IC) and breaking (IB) parts of 
the potential, respectively, and v^j and vfj are the momentum-independent and momentum- 
dependent parts of the interaction. For the AV18, ifij corresponds to the contributions of 
the static components, including isospin-independent and isospin-dependent central, spin- 
spin, and tensor terms, while vfj retains the contributions from the spin-orbit and quadratic 
momentum-dependent components. The v-^^ part in the AV18 is parameterized by the four 
operators 

Oj^j = Tij , <Tj ■ (TjTij , Sij Tij , (tj^j; -|- Tj 2) , (22) 

where S^j is the standard tensor operator and the isotensor operator Tij is defined as Tij = 

The MI two-body currents arising from v^^ have been constructed following the procedure 



of Ref. [l8j, which will be hereafter referred to as the meson-exchange (ME) scheme. Within 
this scheme, the isospin-dependent static part of v^j is assumed to be induced by exchanges of 
effective pseudoscalar (PS), or "vr-like," and vector (V), or"p-like," mesons. The propagators 
associated with these exchanges are projected out of the (isospin-dependent) central, spin- 
spin, and tensor components of v^y The resulting two-body currents satisfy the CCR with 



v^j by construction. Explicit expressions can be found in a number of references (see Ref. [16 
and references therein). 



The currents arising from vfj have been obtained following the procedure of Ref. 19 
which will be referred to as the minimal-substitution (MS) scheme, reviewed and generalized 
m Ref. M. We first note that the isospin operator Tj ■ r,- can be expressed in terms of the 



space-exchange operator {Pij), using the relation 

Ti ■ Tj = -1 - (1 + CTj ■ crj)Pij , (23) 

valid when operating on antisymmetric wave functions. The operator Pij is defined as 
Pij = e^-'''^'"'"^'^'^^ , where the V-operators do not act on the vectors rij = Vi — Vj = —Vji 
in the exponential. In the presence of an electromagnetic field, minimal substitution is 
performed both in the explicit momentum dependence of the two-nucleon potential as well 
as in the implicit momentum dependence implied by Tj ■ Tj. The resulting current operators 
have been derived in Ref. [l6|. Here we only list the final result for the current operators 
associated with the isospin-independent and isospin-dependent spin-orbit interaction. In 
this case, vfj can be expressed as 

4 = + vlijTi ■ Tj , (24) 

where 

vlij = vis{r)L- S , 

V2,ij = visr{r)L- S , (25) 
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with L = Vij X (pi —pj)/2, and S is the total spin of pair ij and vis and visr are the spin-orbit 
parts of the AW potential. Performing minimal substitution in Vi, we obtain 



(26) 



jijiq;ls) = -Vlsir)i^e^e"^^''' - eje^'^'^^j S x r^j 
For the isospin-dependent term ff? ^6 first symmetrize it as 



-Visr{r) {L- STi- Tj +Ti-TjL- S) 



(27) 



and the associated current then reads 
1 



jij{q; 1st) = -Visrir)S x r^j ( r/^e"'-''' - r^ie"^'""^ 



+ ^VisAr)Gliql)iT,XT,)A L S 



I dse"i-^+ I dse"^' L-s\ , (28) 



with rji = G^{qf^) Ti ■ Tj + G^{q'pTi^;^, and ds is the infinitesimal step on the generic path •jij 
(7^j) that goes from position i (j) to position j (i). Since the choice of the two integration 
paths 7ij and Aji is arbitrary, the definition given above for jij{q; 1st) is not unique. How- 
ever, whatever choice is made, the corresponding current will satisfy the CCR with ff by 
construction. The simplest choice for ^jij (7^,) is that of a linear path (LP), which leads to 



(29) 



j^'-iq; 1st) = -VisAr)S x r,, ( r^,e^'^-^' - r/.e"''^ 



+ ^VUr{r)Gl{ql){Ti X Tj), 



L-S) r,,U,{q) + n,U,{q) {L ■ S) 



where fij{q) is defined as 



fijiQ) 



q ■ r 



(30) 

and fij{q = 0) = i. It is interesting to note that, in the limit q 0, the two-body 
current operator derived in the MS scheme becomes path-independent and hence unique [16 



Therefore, for processes involving small momentum transfers, the intrinsic arbitrariness of 
the MS scheme may be of little consequence. 
In earlier works, for example Refs. 



20|, 1211, 122, 123|, 124, the two-body currents from the 



spin-orbit interaction were constructed within the ME scheme, by assuming that its isospin- 
independent components are due to exchanges of "cr-like" and "cj-like" mesons, while the 
isospin-dependent ones originate from "p-like" exchanges. The resulting currents, however, 
are not exactly conserved. This lack of consistency seems to lead to a significant discrepancy 
between theory and experiment in some of the pd radiative capture polarization observables, 
specifically the tensor polarization observables T20 and T21, at low energies [l6|. 

The MI two-body currents arising from the IB terms are generated by the operator 



present in the isotensor operator Tj,-, and are easily constructed [16|. However, their 



contributions have been found to be negligibly small in the study of the electromagnetic 
structure of A=2 and 3 nuclei. This is also the case in the present study of electromagnetic 
transitions in A=Q and 7 nuclei. 

The MD part of the two-body current is purely transverse and therefore is not constrained 
by the CCR. The model adopted here includes the (isoscalar) piT'j and (isovector) uh'J 



6 



TABLE I: Values of the coupling constants Jt^nn-, QpNN-, QujNN-, dfrn-'yi and g^T^-y and monopole form 
factor cutoffs A^^, Ap and A^j used in the present work. 









9piT-y 




A^ 


Ap 


A,^ 


0.075 


0.55 


16.96 


0.56 


0.63 


0.75 GeV 


1.25 GeV 


1.25 GeV 



transition currents, as well as the currents due to excitation of intermediate A isobars. The 
latter are obtained within a nonperturbative treatment based on the transition correlation 
operator approach, reviewed below in Sec. IIIIBi 

The piT'-f and uiT'y MD two-body currents are given by 



I pTT'f 



. fnNNgpNNgpTTj 
1 



X 



■ kj 



(fcj, kj 



. fTTNNgu)NNgujTT-i 



X 



kli 



{kf + ml){B + ml 



{kf + ml){k] + ml 



-T. 



3,z 



(31) 



(32) 



Here fcj (kj) denotes the fractional momentum transfer to nucleon i (j), so that fcj + kj = q. 
The gpn'y, gpNN, gujTT-y, gujNN, and fnNN are the pnj, pNN, ujn'j, uNN, and ttNN coupling 

are the pion, p- and tu-meson masses, respectively. Finally, 



constants. 



while m^r, mp and 



monopole form factors at the pion and vector-meson vertices, given by 



fa{k) 



A2 + F 



n,p,uj 



(33) 



are introduced, to take into account the finite size of nucleons and mesons. The values of 
all the coupling constants and the cutoffs A^ adopted in this work are listed in Table [11 In 
particular, the ottj and uiT'-f coupling constants are obtained from the measured widths of 
the p — s> vr -|- 7 [251] and to ^ ti + j \2&\ decays, while the uNN coupling constant and the 
cutoffs Ajr, Ap and A^ are rather soft but still close to those inferred from models of the NN 
potential. 

In Ref. jl6|, the currents induced by the three-nucleon interaction Vijk associated with 
P-wave two-pion exchange were also constructed. However, their contribution to ^4=3 ob- 
servables was calculated to be quite small. These currents are neglected in the present 
study. 



B. Beyond Nucleons Only 

The simplest description of the nucleus views it as being made up of nucleons, and as- 
sumes that all other sub-nucleonic degrees of freedom may be eliminated in favor of effective 
many-body operators acting on the nucleons' coordinates. The validity of such a descrip- 
tion is based on the success it has achieved in the quantitative prediction of many nuclear 



7 



observables [8| . However, it is interesting to consider corrections to this picture by including 
the degrees of freedom associated with nuclear resonances as additional constituents of the 
nucleus. When treating phenomena which do not involve explicitly meson production, it is 
reasonable to expect that the lowest excitation of the nucleon, the A-isobar, plays a leading 
role. In this approximation, the nuclear wave function is written as 

^iJj^^^ = ^iJ{NN---N) + ¥^\NN---NA) + ¥^\NN---NAA) + -- - , (34) 

where \1/ is the part of the total wave function consisting only of nucleons, '^^^^ is the 
component in which a single nucleon has been converted into a A-isobar, and so on. The 
nuclear two-body interaction is taken as 



B,,B,=N,A b'b'=NA 

J 



(35) 



where transition interactions such as Vij{NN — » NA), Vij{NN — > AA), etc. are responsible 
for generating A-isobar admixtures in the wave function. The long-range part of Vij is due 
to pion-exchange, while its short- and intermediate-range parts, influenced by more complex 
dynamics, are constrained by fitting AW scattering data at lab energy < 400 MeV and 
deuteron properties [27| . 

Once the AW, WA, and AA interactions have been determined, the problem is reduced 
to solving the W -|- A coupled-channel Schrodinger equation. However, this would involve a 
large number of W-|- A channels and therefore the practical implementation of this method is 
very difficult. In a somewhat simpler approach, known as the transition-correlation-operator 
(TCO) method 28|], the nuclear wave function is written as 



N+A 



Kj 



(36) 



where \l/ is the nucleons-only wave function, 5 is a symmetrizer, and the transition operators 



are defined as 



tr 



NA 



AA 



(37) 



NA 
AA 



u 



J ' 

T,;-T 



(38) 
(39) 



Here, and Tj are spin- and isospin-transition operators which convert nucleon i into a 
A-isobar, S-J and 5*/^^^ are tensor operators in which, respectively, the Pauli spin operators 
of either particle i or j, and both particles i and j are replaced by corresponding spin- 
transition operators. The [/*J vanishes in the limit of large interparticle separations, since no 

etc., are obtained 



aril I 



u 



trill 



A components can exist asymptotically. The functions u 
from two-body bound and low-energy scattering state solutions of the full W -|- A coupled 
channel problem, with the Argonne V2% (AV28) model [23] as discussed in Ref. (28| . 

We note that the perturbation theory (PT) description of A-admixtures is equivalent to 
the replacements: 

Vij{NN NA) 



ttNA,PT 
ij 



AA,PT 



m — 



AA) 



2(m — m/\) 



(40) 
(41) 
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where the kinetic energy contributions in the denominators of Eqs. (HOj) and (HT!) have been 
neglected (static A approximation). Note that the transition interactions Vij{NN NA) 
and Vij{NN A A) have the same operator structure as Ujj^ and Ufj^ of Eqs. (!38l) and (!39l) . 
but with the M°"'^"(r) and u^™{r) functions replaced by, respectively, 

(r) = i^^— C(x), (42) 



An 3 X 



-X 



^*™(r) = ^-^^(l + - + ^)'—C\x). (43) 
47r 3 V X X 

Here a = II, III, x = m^r, (//)« = fnNNfnNA, fnNAfnNA for a = II, III, respectively, and 
the cutoff function C{x) = l-e"^^", with A = 4.09. In the AV28 model Q Una = ^Unn- 
This perturbative treatment has been often used in the literature to estimate the effect of A 
degrees of freedom on electroweak observables. However, it may lead to a substantial over 



prediction of their importance |23l . l28l | , since it produces A^A and A A wave functions which 
are too large at short distance. 

The nuclear electromagnetic current is now expanded into a sum of many-body terms as 
in Eq. 0121) . However, here each term operates not only on the nucleon, but also on the A- 
isobar degrees of freedom. Therefore, the one- and two-body currents (ignoring three-body 
currents) are written as 

3l'\q)= E Mq;B^B'), (44) 
b,b'=n,a 

j^iq) = T.s.,b,.n,J:b[,b;.na ^^M^ 1^,1^^ . (45) 

In the present work, however, we only keep the purely nucleonic two-body currents discussed 
in the previous section. 

The one-body N ^ A transition and A currents are given by 

Mq; N^A) = --i-G,^^(gJ)e^'^--q x S,T,,, , (46) 

j,{q; A ^ A) = -^G,AA(g2)e'''-^'q x S,(l + 6,,) , (47) 

where S (0) is the Pauli operator for the A spin 3/2 (isospin 3/2), and the expression for 
ji{q; A N) is obtained from that for ji{q; N ^ A) hj replacing the transition spin and 
isospin operators by their hermitian conjugates. The A^A-transition and A electromagnetic 
form factors, respectively G^na and G-y^A, are parameterized as 

G,NA{ql) = , (48) 



G.AA(gJ) = , "^""^ . (49) 
(1 + qI/^Ia 



Here the A^A-transition magnetic moment /x^tva is taken equal to 3 fijy, as obtained from 



an analysis of 7A^ data in the A-resonance region [29|; this analysis also gives Aata.i = 
0.84 GeV and Ana,2 = 1-2 GeV. The value used for the A magnetic moment /i^AA is 
4.35 fj,N by averaging results of a soft-photon analysis of pion-proton bremsstrahlung data 
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si, 



near the A^"*" resonance 
the nucleon form factor, 
quadrupole transition. 

pion photoproduction amphtude is found to be experimentally small at resonance [31 
neglected is the A convection current. 



and Aaa = 0.84 GeV as in the dipole parameterization of 
In principle, to A excitation can also occur via an electric 
Its contribution, however, has been ignored, since the associated 

Also 



C. Matrix elements 

Matrix elements of the current operator can be written schematically as 

[(^iV+A,/|^^+A,/)(^7V+A,|^^+A,,,)]^/2 ' ^ ) 

where the initial and final state wave functions include both nucleonic and A-isobar degrees 
of freedom. It is convenient to expand \E'iv+A as 

^jv+A = ^ + Ef^S'^ + --- ' (51) 

i<j 

and the matrix element of the current operator becomes 

(^^+A,/ I J I ^^+A,i) = I j(A^only) I ^i) + I I ^i) . (52) 

Here j(A^only) denotes all one- and two-body contributions to j{q) which only involve 
nucleon degrees of freedom, i.e., j(A^only) = j^^'^N ^ N) + j^'^^NN NN). The 
operator j(A) includes terms involving the A-isobar degrees of freedom, associated with the 
explicit A currents j^^\N A), ^'•^•'(A N) and j^^\A A), and with the transition 
operators f/*J. The operator j(A) is illustrated diagrammatically in Fig. [H The terms (a)- 
(g) in Fig. [1] are two-body current operators, while the terms (h)-(j) are to be interpreted as 
renormalization corrections to the "nucleonic" matrix elements (\t'/ | j(A^only) | \E'j), due to 
the presence of A-admixtures in the wave functions. We note that not included in j(A) are 
all remaining connected three-body contributions of the type of Fig. [2|, which are expected 
to be significantly smaller than those considered in Fig. [TJ 

The terms in Fig. [T] are expanded as operators acting on the nucleons' coordinates. For 
example, the terms (a) and (e) have the structure, respectively, 

(a) = ji'\A^N)U$\ (53) 
(e) = t/^t'^Sf (A-A)f//^^, (54) 



which can be reduced to operators involving only Pauli spin and isospin matrices by using 
the identities 

S^-AS-B = ^A-B-^o--(AxB) , (55) 

S^-AS-BS-C = -iA-(BxC)-icr-AB-C 

3 3 

-iA-BC-cr + ^A- (B-cr)C , (56) 
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FIG. 1: Diagrammatic representation of operators included in j(A) due to one-body currents 

jW{N A), j(i)(A ^ iV) and j(^)(A ^ A), and transition correlations U^^, U^^ , U^^, and 

corresponding hermitian conjugates. Wavy, thin, thick, dashed and dashed with a x lines denote 

' 't 

photons, nucleons, A-isobars, and transition correlations U^^ and U^^ , respectively. 



(a) 




(b) 





(c) 



(d) 



(e) 



(f) 



(g) 



(h) 



(i) 



a) 



FIG. 2: Diagrams associated with connected three-body terms, which are neglected in the present 
work. Wavy, thin, thick, dotted, and dashed and dashed with a x lines denote photons, nucleons, 

A-isobars, the two-body current j^'^\NN — > NN), and the transition correlations U^^ and U^^ 
respectively. 



'X- 



(a) 



(b) 



i» • • 



(c) 



where A, B and C are vector operators that commute with cr, but not necessarily among 
themselves. Expressions for the other terms of Fig. [T] are obtained in a similar fashion. 

The denominator of Eq. ( l50l) requires the calculation of the initial and final state wave 
function renormalizations, which are given by 



A\2 



(iV-) 



N+A 



VP) 



+ 



three — body terms) 



(57) 



and the three-body terms have been neglected consistently with the approximation intro- 
duced in Eq. f l52|) . as discussed above. The wave-function renormalizations {N^) for the 
different nuclei considered in the present work are listed in Table [TTl The TCO approxi- 
mation of Eq. (|5T|) gives a renormalization of 1.0023 for the deuteron, whereas the exact 
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TABLE II: Wave function renormalizations, A^'^, obtained for the A=2-7 nuclei, when the TCO 
calculation is based on the AV28 interaction is used with purely nucleonic GFMC wave functions 
for the AV18+IL2 Hamiltonian model. 



2H 


3H 


^He 


6Li 


^Li 


■^Be 


1.0023 


1.016 


1.016 


1.050 


1.071 


1.073(1) 



coupled-channel result for the AV28 potential is 1.0026. Note that PT estimates of the im- 
portance of A-isobar degrees of freedom in photo- and electro-nuclear observables typically 
include only the contribution from single ^ A transitions (namely diagrams (a) and (b) 
in Fig. [1]) and ignore the change in the wave function normalization. 



IV. RESULTS 

In Ref. [l| we reported results for fifteen electroweak transitions in A=6,7 nuclei. We 
pointed out that MEC contributions are expected to be significant in magnetic moment 
(/i) and magnetic dipole (Ml) transition calculations. Here we present our results for fi for 
A=2-7 and Ml transitions for A=6,7 including the MEC contributions. The first subsection 
discusses the magnetic moment results for A=2-7 nuclei and the second subsection discusses 
the Ml transitions in A=6,7 nuclei. 

The calculation of /i's or Ml transitions is fairly straightforward. For example, fi is 
obtained from the diagonal matrix element 

2 TTl 

fi = -i\im—{r,Mj = J;T \ jyiqx.) \ J\Mj = J;T) , (58) 

where the momentum transfer q is taken along the x axis, the nuclear state with angular 
momentum quantized along the z axis has Mj=J, and m is the nucleon mass. The VMC or 
GFMC wave function for the given ( J'^, Mj = J; T) state is then constructed. Evaluation 
of the various contributions is done as a function of the momentum transfer q for several 
small values q < 0.05 fm"^ and then extrapolated smoothly to the limit q=0. The error due 
to extrapolation is much smaller than the statistical error in the Monte Carlo sampling. 

In the tables below, we present the one-body, i.e. impulse approximation (lA), results and 
contributions from various pieces of the two-body MEC operators separately: pseudoscalar + 
vector (PS+V), minimal-substitution (MS), model-dependent (MD), and A. The A column 
includes contributions from both the explicit MEC terms of Eqs. (H6l) and fH7|) and the 
renormalization of the nucleons-only terms as given by Eq. (!50l) : 

^ ^ {^f\jiA)+jiN only) I VI/,) _ {^fljjN only)|vl/,) 



A. Magnetic Moments in A=2—7 Nuclei 

Table IIIII shows the magnetic moment results for ^4=2-7 nuclei. The last two columns 
list the total magnetic moments and corresponding experimental numbers [s^ . 33l |. The 
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TABLE III: Magnetic moments in nuclear magnetons for A =2-7 nuclei; lA, PS, V, MS, MD stand 
for impulse approximation, pseudoscalar, vector, minimal-substitution, and model-dependent, re- 
spectively. Details can be found in the text. IS and IV in the "Nucleus" column denote the 
isoscalar and isovector combinations of the preceding nuclei. The experimental values are from the 
compendia [s^] except for the very recent measurement for ^Be {ssl; they have been rounded to 3 
decimal digits, except for ^H. 



Nucleus 


Method 


lA 


r kj-f- V 


MEG 

iViO 




A 
LA 


Total 


Expt. 






0.8467 





-0.0022 


0.0031 


0.0009 


0.8485 


0.8574 


n 


V iVlU 


Z.OoU 


n Q 1 Q 


n nn9 
-u.uuz 


n ni 7 

U.Ui ( 


n ni H 
U.Uio 


z.yoz^i ) 


9 Q7Q 

z.y (y 


1 1 




Z.O ( O^ZJ 


yj.oZZyz, ) 


n nn9 
-u.uuz 


n ni 7 

U.Ui ( 


n ni /I 

U.Ui'i 


9 QOAt'i\ 

z.yz^^^o ) 


9 Q7Q 

Z.y / y 


11 


ri rl 




n Q91 


-U.UUi 


n ni 7 
U.Ui ( 


n ni /I 
U.Ui4 


z.yzo 


9 n7n 
z.y i\) 


iie 


V iviv^ 


-i . / DO 


-U.Oi ( 


n nm 

-U.UUi 


n ni n 

-U.UiU 


n ni 

-U.Uio 


9 1 r\7( 1 

-Z.iU 1 


9 1 98 
-Z. iZo 


rle 




1 v^^^/'9^ 
-i. ( 0\)[Z ) 


-yj.oLoyz ) 


n nni 
-U.UUi 


n ni n 
-U.UiU 


n ni Q 
-U.Uio 


9 1 nQ f Q\ 

-Z.Wo[o ) 


9 1 9C 

-Z. iZo 


^He 


HH 


-1.764 


-0.316 


-0.001 


-0.010 


-0.014 


-2.105 


-2.128 


IS 


GFMC 


0.408 


0.001 


0.002 


0.003 


-0.003 


0.411 


0.426 


IV 


GFMC 


4.329 


0.640 


0.001 


0.027 


0.030 


5.027 


5.107 


6Li 


VMC 


0.815(1) 





-0.008 


0.004 


-0.006 


0.805(1) 


0.822 


6Li 


GFMC 


0.810(1) 





-0.007 


0.005 


-0.008 


0.800(1) 


0.822 


^Li 


VMC 


2.906(4) 


0.318(3) 


-0.011(1) 


0.019 


-0.042 


3.190(7) 


3.256 


^Li 


GFMC 


2.870(8) 


0.340(6) 


-0.009(4) 


0.020 


-0.053 


3.168(13) 


3.256 


■^Be 


VMC 


-1.098(5) 


-0.317(6) 


0.005 


-0.012 


-0.078 


-1.503(5) 


-1.400 


^Be 


GFMC 


-1.058(9) 


-0.343(6) 


0.007 


-0.011 


-0.088 


-1.493(15) 


-1.400 


IS 


GFMC 


0.906 


0.001 


-0.002 


0.004 


-0.073 


0.836 


0.929 


IV 


GFMC 


3.928 


0.683 


-0.019 


0.030 


0.039 


4.661 


4.654 



total fi is obtained from the sum of the lA and the two-body contributions from various 
pieces. Apart from the case, we present the VMC results followed by the GFMC results 
in the following row for each magnetic moment. Hyperspherical harmonics (HH) results 
for the trinucleons interacting by the AV18 AW potential and older Urbana-IX (UIX) [s^] 
3N potential are shown for comparison. We also present the GFMC isoscalar and isovector 
combinations for ^4=3, 7. 

The results presented in Table IIIII show the significant impact of the two-body operators 
in those cases with nuclear isospin T = 1/2. MEC contributions boost the lA by about 
16% in the ^4=3 isovector case and by about 19% in the A=7 nuclear states. For the two 
T = states, namely ^H and ^Li, we see that the lA magnetic moments are not modified 
significantly by the MEC, as expected for any isoscalar state. We note, however, that the 
present isoscalar MEC contributions to the deuteron magnetic moment are smaller than 
reported previously in Ref. [2|. This is the result of the different way in which two-body 
currents from the momentum dependent components of the AV18 have been constructed in 
this work (see Ref. [l6| for a discussion of this issue). 
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The magnetic moments for and ^He are closer to the experimental values of /i in both 
cases when the MEC contributions are added to the lA values. The major contributions 
come from the pseudoscalar and vector piece of the currents. The model-dependent piece 
and A contributions are small, but not insignificant, while the minimal substitution piece is 
tiny. We also note that the VMC, GFMC, and HH results are all very close to each other 
for all the separate pieces as well as for the total /i. 

The simplest picture of consists of a S" = pair of neutrons and a proton, all in a total 
L = state. That of ^He is the same with proton and neutron interchanged. In this picture 
the impulse approximation magnetic moments of and '^He would be the same as those 
of the proton and neutron respectively: 2.79 and -1.91. As can be seen in Table llllt the 
GFMC impulse magnetic moments are 8% smaller in magnitude than these. This is largely 
due to the tensor force which has several effects: 1) the wave functions contain ~ 9% of 
L = 2 components which result in orbital contributions of +0.04 and +0.05 to the magnetic 
moments of and ^He, respectively; 2) the odd nucleon is not 100% aligned with the 



nuclear spin; and 3) the pair of like nucleons has ~ 10% 5* = 1 component [35j. The last 
two effects reduce the spin contributions to the magnetic moments from the pure nucleon 
values to 2.53 and -1.81. 

A similar analysis of the A=7 magnetic moments can be made. The ground state of 
^Li looks a great deal like an a particle plus triton in relative P-wave motion, whereas 
''Be looks like a plus ^He. Thus, the orbital contribution to the A=7 impulse magnetic 
moments is expected to be significantly larger than for A=3; it is 0.42 and 0.67 for ^Li 
and ^Be, respectively. The spin contributions, 2.43 and -1.72 are quenched relative to the 
spin contributions to the A=3 values; this is probably due to tensor interactions between 
a nucleon in the a core and one in the valence A=3 cluster. The MEC contributions are 
generally larger in the A=7 nuclei, again probably due to interactions between the core and 
valence nucleons. 

The GFMC propagation for the magnetic moment of the ^Li(| ) ground state is shown 
in Fig. [3l In this figure the two solid purple lines correspond to the VMC (impulse and 
total) values, the green squares (red circles) are extrapolated GFMC impulse (total) propa- 
gations and the solid green (red) lines starting at r=0.1 MeV~^ are the final GFMC averages 
with dashed lines to denote the Monte Carlo error. The GFMC propagation reduces the 
VMC matrix element slightly in both cases. The GFMC impulse result is only 88% of the 
experimental value, but the MEC contributions raise the total to 97%. 



B. Magnetic Dipole Transitions in A=6,7 Nuclei 

In Table HV] we present the different contributions to the matrix elements for Ml tran- 
sitions in A=Q,7 nuclei. As in the case of magnetic moment calculations, we see that the 
most significant contributions come from the pseudoscalar and vector pieces of the two-body 
current operators. 

The first two rows of Table IIVI show the various pieces of the Ml matrix element for 
^Li(0+;1) -^^Li(l+;0) transition. We note that both VMC and GFMC lA results are 
boosted by 7-8% by MEC. The corresponding decay widths for this transition are pre- 
sented in Table El The total width we obtain from the VMC calculation agrees very well 
with the experimental value, whereas the total GFMC width is slightly outside of the present 
experimental range. Figure H] shows the matrix elements for this function of r. 

As in the previous figure, the two solid purple lines represent the VMC impulse and total 
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FIG. 3: (Color online) Extrapolated GFMC magnetic moment for the '^Li(| ) ground state in 
impulse approximation (green squares) and with MEG (red circles). VMC values (purple lines) 
and averaged GFMC values (lines with error bars) are also shown. 



estimates, the green squares represent the GFMC propagated points for impulse and the 
red circles denote the total GFMC matrix elements. The average GFMC results are shown 
as solid lines with error bars starting at r=0.1 MeV^^. All the GFMC points as well as the 
averages represent the extrapolated matrix elements. We note that, also in the present case, 
the GFMC propagation slightly decreases the VMC value for the ^Li(0+; 1) ^6Li(l+;0) 
transition. We see that the propagated points are quite stable with r, which we ran up to 
0.3 MeV~^ in this case. 

Table IVl also shows two Ml transitions in A=7 nuclei. The MEG corrections are 15-17% 
of the Ml matrix elements obtained in both VMC and GFMC calculations. The model 
independent pieces (PS+V) are the largest contributions and are also very similar for both 
^Li and ^Be. The decay widths for ^Li(i") ^'^Lid") and ^Be(l") ->'^Be(|") Ml transitions 
are shown in Table El The total widths from both VMC and GFMC match very well with 
the experimental decay widths. 
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TABLE IV: Matrix elements of Ml transitions in A=6,7 nuclei. Column headings are defined in 
Table imi 



TIT , TTT 

Jj Jf 




Method 


lA 




MFC 




Total 








PS+V 


MS 


MD 


A 




6Li(0+;l) 


^6Li(l+;0) 


VMC 


3.683(14) 


0.307 


0.003 


0.010 


-0.053 


3.950(14) 


6Li(0+;l) 


^6Li(l+;0) 


GFMC 


3.587(16) 


0.323 


0.002 


0.012 


-0.048 


3.876(14) 


^Li(r) - 


.^Li(i-) 


VMC 


2.743(17) 


0.396 


0.006 


-0.017 


-0.034 


3.162(22) 


^Li(i-) - 


.^Li(i-) 


GFMC 


2.677(19) 


0.395 


0.011 


-0.017 


0.072 


3.138(22) 


'Be(r) - 


-^Be(i-) 


VMC 


2.420(30) 


0.390 


-0.005 


0.010 


-0.024 


2.791(36) 


'Be(i-) - 


-^Be(|-) 


GFMC 


2.374(31) 


0.394 


-0.010 


0.010 


-0.002 


2.766(36) 


TABLE V: Impulse approximation (lA) and total Ml transition widths in 


eV for A= 


=6,7 nuclei. 


TTT , TTT 

Ji Jf 




Mode 


VMC 




GFMC 


Expt. 






lA 


Total 


lA 


Total 




6Li(0+;l) 


^6Li(l+;0) 


Ml 


7.09(6) 


8.15(6) 


6.72(6) 


7.85(6) 


8.19(17) 


^Li(i-) - 


.^Li(i-) 


Ml (10- 


-3) 4.75(6) 


6.31(9) 


4.52(6) 


6.21(9) 


6.30(31) 


'Be(i-) - 


-^Be(r) 


Ml (10- 


-3) 2.62(7) 


3.49(9) 


2.52(7) 


3.42(9) 


3.43(45) 



V. CONCLUSIONS 

In summary, we have reported results for the magnetic moments and magnetic dipole 
transitions in nuclei with mass numbers A < 7. The calculations have used essentially 
exact wave functions derived from a realistic Hamiltonian that reproduces well the low-lying 
spectra of these nuclei as well as of those in the mass range ^4=8-10. Leading terms in 
the nuclear electromagnetic current have been constructed to satisfy current conservation 
with the two-nucleon potential, AV18, used in the Hamiltonian. Additional contributions 
associated with the explicit presence of A isobar degrees of freedom have been accounted 
for by including, in an approximate fashion, A components in the nuclear wave functions 
with the transition-correlation-operator method. 

Overall, the agreement between the calculated and experimental magnetic moments and 
transition rates is quite satisfactory, particularly in the isovector channel where differences 
between computed and experimental amplitudes are ~ 1.5%. On the other hand, in the 
isoscalar channel these differences seem to progressively become worse as the mass number 
increases; they are about 1% in deuteron, 4% in ^He/^H, and 10% in ''Be/^Li. Of course, 
isoscalar transitions are suppressed both at the one- and two-body levels: the lA current is 
proportional to the nucleon isoscalar magnetic moment, which is five times smaller than the 
corresponding isovector combination; leading two-body currents from pion-exchange and A 
excitation have isovector character. Two-body isoscalar contributions arise in the present 
study from short-range mechanisms: the momentum-dependent components of the AVIS, 
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X (MeV ^) 



0.30 



FIG. 4: (Color online) Extrapolated GFMC Ml matrix element for the '^Li(0+;1) ^'^Li(l+;0) 
transition in impulse approximation (green squares) and with MEC (red circles). VMC values 
(purple lines) and averaged GFMC values (lines with error bars) are also shown. 

the p7r7 transition current, and renormalization corrections induced by A admixtures in the 
wave functions, Eq. fl59|) . 

We conclude by noting that in a chiral effective-field-theory framework isoscalar correc- 
tions are suppressed by {Q/A^Y {Q denotes a generic small momentum and A^^ ~ 1 GeV is 



36, 37 



the chiral-symmetry-breaking scale) relative to the leading-order (LO) lA current 
These N^LO corrections have been calculated in the deuteron and trinucleon isoscalar mag- 
netic moments, and are ~ 1% relative to LO but of opposite sign, so that they increase the 
discrepancy between theory and experiment. At N^LO, or (Q/A^^,)^, a number of isoscalar 
two-body currents originate from four-nucleon contact interactions involving two gradients 
of the nucleon fields ^]. Their contributions to electromagnetic observables have yet to 
be calculated. It will be interesting to see whether these isoscalar currents as well as the 
corresponding isovector ones up to N'^LO will improve the present picture. 
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